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ABSTRACT 


We investigate the effect of magnetic reconnection between open and closed field, 
(often referred to as “interchange” reconnection), on the dynamics and topology of 
coronal hole boundaries. The most important and most prevalent 3D topology of the 
interchange process is that of a small-scale bipolar magnetic field interacting with a 
large-scale background field. We determine the evolution of such a magnetic topology by 
numerical solution of the fully 3D MHD equations in spherical coordinates. First, we 
calculate the evolution of a small-scale bipole that initially is completely inside an open 
field region and then is driven across a coronal hole boundary by photospheric motions. 
Next the reverse situation is calculated in which the bipole is initially inside the closed 
region and driven toward the coronal hole boundary. In both cases we find that the stress 
imparted by the photospheric motions results in deformation of the separatrix surface 
between the closed field of the bipole and the background field, leading to rapid current 
sheet formation and to efficient reconnection. When the bipole is inside the open field 
region, the reconnection is of the interchange type in that it exchanges open and closed 
field. We examine, in detail, the topology of the field as the bipole moves across the 
coronal hole boundary, and find that the field remains well-connected throughout this 
process. Our results imply that open flux cannot penetrate deeply into the closed field 
region below a helmet streamer and, hence, support the quasi-steady models in which 
open and closed flux remain topologically distinct. Our results also support the 
uniqueness hypothesis for open field regions as postulated by Antiochos et al. We discuss 
the implications of this work for coronal observations. 


Subject Headings: Sun: corona - Sun: magnetic fields - Sun: reconnection - Sun: coronal 
hole 
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1. Introduction 

The solar magnetic field is the primary driver of solar activity and is the principal 
conduit for coupling the energy of the sun’s convective envelope to the corona and, 
subsequently, to the solar wind. The question of the structure and dynamics of the coronal 
magnetic field is, therefore, central to understanding all solar activity and the structure 
and evolution of the heliosphere. With the advent of XUV/X-ray imaging from space 
missions, it became possible to observe coronal structure directly. Even the low- 
resolution images from the early SKYLAB mission showed clearly that the large-scale 
corona is composed of two physically distinct regions: “closed-field” regions, consisting 
primarily of bright X-ray loops, and “coronal holes” that are dark in X-rays (Zirker, 

1977). The photospheric flux below coronal loops is observed to be bipolar implying that 
the field is closed, i.e., connected to the photosphere at both ends of the coronal field 
lines. On the other hand, the photospheric flux below coronal holes is unipolar, on 
average, implying that the field there is open, i.e., connected to the photosphere at one 
end with the other extending outward indefinitely into the heliosphere. Coronal holes, 
therefore, are a source region for solar wind, which also explains why these regions are 
dark in X-rays. The coronal density is low there due to the large mass and energy flux 
required to power the solar wind. 

Motivated by the XUV/X-ray observations and by the basic theory of the solar 
wind given by Parker (1958), a standard model has developed for the large-scale coronal 
magnetic field, the “quasi-steady model” (e.g., Antiochos et al. 2007). The underlying 
assumptions of this model are that the magnetic field is determined by the instantaneous 
distribution of the large-scale radial flux at the photosphere and the balance between gas 
pressure and magnetic stress in the corona. In the quasi-steady model, the coronal 
magnetic field is assumed to be static with a smooth structure consisting of topologically 
well-separated open and closed field regions. Of course, the flux distribution at the 
photosphere, even the large-scale distribution, does change with time due to flux 
emergence and photospheric motions; but, since typical Alfven speeds over much of the 
corona are orders of magnitude greater than the driving photospheric flows, the coronal 
field evolution can be represented as a series of time stationary states. Note that implicit 
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in the quasi-steady model is the assumption that flux opens and closes, (most likely 
involving reconnection), in response to the photospheric changes. 

The first and simplest implementation of the quasi-steady model was the potential 
field source surface model (Altschuler & Newkirk 1969; Schatten et al. 1969; Hoeksema 
1991). In this model gas pressure is neglected and the magnetic field is taken to be 
current-free. The open flux is determined by the assumption that the field is purely radial 
at some given radius (the source surface). Although the assumptions of the source surface 
model are extreme, it has proved to be highly useful, because it is easy to calculate and is 
surprisingly accurate at reproducing the observed pattern of open and closed regions on 
the Sun (e.g., Hoeksema 1991). Over the last decade, numerical models have been 
developed that calculate steady-state solution to the full magnetohydrodyanmic (MHD) 
equations, so that the assumptions in the source surface model can be relaxed (e.g., 

Linker et al. 1999; Odstrcil 2003; Roussev et al. 2003). In fact, these models can even 
drop the quasi-steady assumption and include the photospheric time dependence, but this 
is rarely done due to the added complexity and computational expense. Furthermore, a 
rigorous time-dependent model would require robust treatment of flux emergence and 
cancellation, which is not yet available. 

In addition to capturing the observed distribution of coronal holes on the Sun, the 
quasi- steady models are fairly accurate in reproducing in situ measurements of the steady 
solar wind magnetic field and plasma (e.g., Zurbuchen 2007; Lepri et al. 2008). 
Furthermore, the dynamics implicit in the model are in qualitative agreement with 
coronal plasma observations. The observation of plasma inflows and outflows, (e.g., 
Hundhausen et. al., 1984; Howard et. al., 1985; Sheeley & Wang, 2002), the observation 
of quasi-rigid rotation of coronal holes, and the existence of the highly variable slow 
wind suggest continuous opening and closing down of flux at the heliospheric current 
sheet, as predicted by the model. 

There is one heliospheric observation, however, that appears to be in direct 
conflict with the quasi-steady model - the measurement of electron heat flux in the solar 
wind. In order to close down heliospheric flux, reconnection between two open field lines 
must occur at an altitude below the Alfven point, where the magnetic energy still exceeds 
the thermal energy. Such a reconnection will create two loops: one having both foot 
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points anchored to the solar surface remaining below the Alfven point, and the other - an 
inverted-loop - entirely detached from the Sun and dragged away with the solar wind. It 
is exactly this type of reconnection process that is implied by coronal observations of the 
streamer belt evolution. Conversely, the opening of previously closed flux requires that a 
loop expand into the heliosphere and be dragged outward by the wind. It has long been 
recognized that such processes should produce a signature in the field-aligned 
suprathermal electron beams (~70 eV to several keV) in the heliosphere (Gosling 1990). 
Streaming electron beams directed away from the hot corona, indicate open flux attached 
at a single foot point. Field lines with both foot points anchored in the solar surface and 
dragged into the heliosphere by the solar wind would exhibit bi-directional, counter- 
streaming electrons, whereas inverted-loops would be devoid of these suprathermal 
electron beams altogether, a so-called “heat flux dropout”. Thus, the heat-flux electrons 
provide a local measure of the global field-line field topology, and as such, are a 
predictive indicator of flux opening and closing. The key inconsistency between 
heliospheric observations and the quasi-steady model is that bi-directional electron beams 
and heat-flux dropouts in the solar wind are rarely observed outside interplanetary 
coronal mass ejections (CMEs) (McComas et al. 1989, 1991; Lin & Kahler, 1992; Pagel 
et al. 2005). 

Motivated by these electron observations, which imply negligible field line 
opening or closing, (and by in situ observations implying large field line wandering in the 
heliosphere), Fisk and co-workers have proposed an alternative theory for the 
solar/heliospheric magnetic field: the “interchange model” (Fisk et al. 1999; Fisk & 
Schwadron 2001; Fisk 2005; Fisk & Zurbuchen 2006). In this model the open flux is 
assumed to be constant during a solar cycle, except for the transient flux of CMEs. This 
assumption appears to be well supported by observations, which show only small 
variation from cycle to cycle in the total heliospheric flux at solar minimum when the 
effect of CMEs can be accurately removed from the heliospheric data. Note, however, 
that the observations for the latest minimum, cycle 23, lower the minimum heliospheric 
open flux threshold, and may even contradict the constant-flux assumption altogether 
(Fisk & Zhao 2009). 
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The dominant process determining the open-field’s evolution in the interchange model is 
magnetic reconnection between open and closed flux, which always conserves the 
amount of each type of flux (e.g., Crooker et al. 2002). It should be emphasized that the 
reconnection postulated by the interchange model is quite different than that in the quasi- 
steady model. In the latter, reconnection occurs primarily at the heliospheric current 
sheet, because that is where open field lines close down. Although field line opening does 
not require reconnection, the opening often involves the ejection of a plasmoid from the 
top of a streamer, which implies reconnection again at a newly-formed heliospheric 
current sheet. Reconnection in the interchange model, on the other hand is statistical in 
nature, and occurs primarily between open flux and the closed flux of coronal loops 
leading to a diffusive motion of open field in the low corona. The open field is derived to 
mix indiscriminately with the closed, throughout the corona, so that reconnection 
between open and coronal-loop field occurs continuously. The interchange model, 
therefore, postulates a very different magnetic topology than the well-separated open and 
closed topology of the quasi-steady. We conclude, therefore, that these topological 
differences are a strong discriminator between the two models of the open field. For the 
interchange model to be valid the open-field topology must be discontinuous and 
inherently dynamic, whereas for the quasi-steady to hold, the topology must remain 
continuous throughout any coronal-field evolution. 

In recent work we analyzed the topological properties of the quasi-steady models 
and derived severe constraints on the possible structure of the open field. We derived, in 
particular, the uniqueness conjecture, which states that irrespective of the complexity of 
the photospheric flux distribution, every unipolar region on the photosphere can contain 
at most one coronal hole (Antiochos et al. 2007). Note that such a topology in which the 
open field has well-defined, connected structure is the exact opposite of that of the 
interchange model. 

The validity of the uniqueness conjecture and of the quasi-steady model, in 
general, turns out to depend critically on the properties of interchange reconnection. The 
key point is that reconnection between open and closed flux is expected to be a generic 
feature of the solar corona and, therefore, must be incorporated into all coronal models, 
including the quasi-steady. Due to the so-called magnetic carpet (Schrijver et al. 1997), 
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coronal holes are obviously not unipolar; they contain numerous small bipoles and, 
therefore, closed flux. As these bipoles move with the photospheric flows, they will 
interact with the open field and undergo interchange-type reconnection. In order for the 
quasi-steady assumption to remain valid during such open-closed interactions, the 
magnetic topology must remain smooth, with the open and closed flux topologically well 
separated. Reconnection, however, requires the formation of current sheets, which are 
topological discontinuities, and reconnection generally gives rise to strong dynamics, 
which would invalidate the quasi-steady assumption. Consequently, it is not clear that the 
magnetic topology would remain smooth during actual time-dependent interchange 
reconnection. Our first objective in this paper, therefore, is to calculate the rigorous 3D 
evolution of a closed field bipole as it moves through and interacts with open field and 
determine whether the resulting structure and dynamics are compatible with the quasi- 
steady assumptions or whether the topology becomes discontinuous as in the interchange 
model. 

A related and equally important issue is the interaction of the closed field of a 
bipole with a coronal hole boundary. In our analysis of the quasi-steady model, we found 
that this type of interaction plays a central role in determining the coronal topology, 
including uniqueness and several other properties (Antiochos et al. 2007). We argued in 
that work that reconnection would enforce the uniqueness constraint, but this was only a 
conjecture. The second objective of this paper is to calculate the time-dependent 
dynamics of coronal hole boundaries rigorously and test the conjectures. We describe 
below our numerical simulations of 3D interchange reconnection and the conclusions for 
coronal structure and dynamics. It should be emphasized that since the physical systems 
we calculate are very general and expected to be ubiquitous on the Sun, our results are 
important for understanding not only the quasi-steady, but any model for the coronal 
magnetic field, including the interchange. 

2. The Topology of 3D Interchange Reconnection 

In order to perform a rigorous study of interchange reconnection, we first need a 
physically robust magnetic field. This requires a quantitative description of a complete 
field topology, not simply a largely two-dimensional sketch of a few open and closed 


-7 - 



field-lines as used in many previous studies. The simplest and most common magnetic 
configuration that can describe interchange reconnection is that of a global bipolar field 
with open and closed regions, and a small-scale closed bipolar region. We can calculate 
this field exactly with an analytic source-surface model that uses the method of images 
(Antiochos et al. 2007). The scalar potential for the source surface field due to a global 
dipole at Sun center and an arbitrary number of smaller dipoles below the solar surface is 
given by: 



where Mi and Tj are the magnetic dipole and position vectors respectively of dipole 
source i, and R is the source surface radius. Note that, for simplicity, we have taken the 
global dipole M 0 at Sun center (x = 0) to be vertical, parallel to the polar axis, and the 
smaller dipoles M; to be horizontal, perpendicular to the radius vector. Their orientation 
in the horizontal plane, however, can be arbitrary. From this potential, the magnetic field 
in the volume is obtained directly from B = grad <E>, and as can be verified by 
straightforward calculation, is purely radial at the source surface, r = R. 

Although the formula above can be used to describe fields of arbitrary 
complexity, the fundamental topology of interchange reconnection is most clearly seen 
by focusing on the case of a global dipole and a single near-surface dipole. Such a field 
is shown in Figure 1, for the source surface position at R = 3 Rs, a Sun-center dipole of 
strength IMol = 10 G oriented toward polar north. The active-region bipole is positioned 
at Iril = 0.9 Rs (below the photosphere) and 49.5° latitude (i.e., north of the equator inside 
the northern coronal hole), with a magnitude IM;I = 50 G oriented along the surface (i.e., 
with no radial component) toward the south pole. As expected, the global dipole produces 
a large-scale, axi-symmetric coronal magnetic field consisting of polar coronal holes and 
closed flux at lower latitudes (Figure 1). The near-surface dipole produces a small bipolar 
flux distribution that, for the particular parameters selected, is completely inside the 
northern, positive-polarity coronal hole. Figure lb shows a close-up of the photospheric 
flux distribution in the hole. Note the presence of the closed polarity inversion line 
surrounding the negative-polarity region of the bipole. The field near this polarity 
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inversion line is low-lying and must close across it; consequently, there must be some 
closed flux inside the coronal hole. This is true irrespective of the size of the negative 
polarity region. There must be a closed field region associated with every bipole in a 
coronal hole. 

Figure lb shows the coronal magnetic field above the small bipole. Its structure 
consists of a hemispherical volume of closed flux surrounded by a background of open 
coronal-hole flux. The closed-field is topologically separated from the open by a dome- 
shaped surface. This topology is simply that of the well-known embedded bipole with its 
fan surface, spine lines, and null point (e.g., Greene 1988; Lau & Finn 1990; Antiochos 
1990; Priest & Titov 1996). The intersection of the fan surface with the photosphere 
forms a closed separatrix curve that defines the boundary between the flux that closes 
across the polarity inversion line to that connecting to the source surface. In other words, 
this photospheric separatrix curve is a coronal hole boundary. All the field lines whose 
photospheric footpoints lie on this curve can be considered to converge onto the null 
point, where they split into the inner and outer spine lines. It should be emphasized that 
although the topology is discontinuous at the fan and spines (i.e., the magnetic 
connectivity is clearly multi-valued there), the magnetic field itself is smooth 
everywhere. In fact, the formula above yields a potential field that is analytic everywhere 
in the interior of the volume. Furthermore, there is no mixing of open and closed field. 

All of the flux inside the fan is closed, whereas all of the flux outside is open. The fan, 
itself, is a singular surface as with every coronal hole boundary in that the fan field lines 
split at the null, so they can be considered to both be open and closed. 

The field of Figure 1 is the fundamental topology in which interchange 
reconnection takes place. It is, by far, the most common multi-polar magnetic topology 
on the Sun, because it is present whenever a parasitic polarity region on the photosphere 
occurs inside some larger, unipolar flux. This topology is expected for essentially every 
magnetic carpet, or larger, bipole on the photosphere. Numerous observations show clear 
evidence for this topology in coronal holes (e.g., Golub et al. 1974), and both potential 
and force-free extrapolations of almost every observed photospheric flux distribution find 
this topology in both open and closed magnetic regions (Aulanier et al. 2000; Fletcher et 
al. 2001; Luhmann et al. 2003). 
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The key question for the coronal models is whether the embedded bipole topology 
remains smooth, with well- separated regions, once photospheric motions stress the field 
so that closed and open lines interact via interchange reconnection. A rigorous answer to 
this question requires solution of the fully-dynamic equations, as presented below, but we 
claim that considerable insight can be obtained by considering the heuristic model for the 
stressing and reconnection illustrated by Figure 2. There are two basic assumptions 
underlying this model. First, we can separate the ideal and resistive response of the 
system so that it evolves, first, purely ideally to some quasi-equilibrium, and then it 
relaxes by reconnection. This approach is not without justification, because reconnection 
will not begin until the system has formed substantial current sheets. The second 
assumption is that the two flux systems on either side of the fan surface move 
independently of each other, except that they always share a common boundary, the fan 
surface, which itself is free to deform. Again this assumption has justification; since the 
photospheric connectivity is discontinuous at the fan, the magnetic stresses due to 
photospheric driving will be discontinuous there, which will give rise to discontinuous 
motions in the corona. Note that even if viscosity were included in the system so that no 
true discontinuity forms, we would expect the gradients of the motions across the fan to 
grow exponentially in time and, consequently, the currents there to reach the dissipation 
scale rapidly. 

We can use this model to determine how the embedded-bipole topology would 
respond to a simple footpoint motion that displaces the closed flux system bodily to the 
right, while keeping the open flux more-or-less fixed (Figure 2b). For such a stressing, 
we expect that, during the ideal response, the inner spine line connecting to the parasitic 
polarity dislocates from the outer spine connecting to the source surface. Since each spine 
line “fans out” at the null to form its own surface, such a dislocation implies that the fan 
surface separates into two surfaces that are in contact everywhere, but with field lines that 
are misaligned. The effect of dislocating the spine lines and fan surfaces, therefore, is to 
deform the null point into a 3D null patch and to form a 3D current sheet at the fan. If the 
system were purely ideal then, in principle, it could achieve an equilibrium state 
containing these discontinuities. 
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A small resistivity can now be included to consider the subsequent evolution due 
to reconnection. The system will attempt to relax, as much as possible, back to the 
potential state to minimize its energy. In particular, reconnection at the null-patch can 
destroy the current sheets and, as illustrated in Figure 2c, deform the null patch back to a 
point, thereby realigning the spines, and if possible, the fans. Note that the evolution just 
described is nothing more than the 3D generalization of Syrovatskii’s (1981) classic 
current sheet formation and null-point reconnection theory (e.g., Antiochos 1996). 

The arguments above suggest that the topology resulting from reconnection will 
maintain clearly separated open and closed field, as in the initial state. A key point, 
however, is that since reconnection conserves any helicity injected to the system by the 
photospheric motions, it cannot undo the photospheric motions and bring the system back 
to a purely potential field. In the evolution illustrated by Figure 2c, the spine lines do not 
actually move, instead different flux tubes become the spines because of interchange 
reconnection. It may be, therefore, that the lowest energy state available to the system 
under the helicity constraints is one with long-lived (up to a dissipation time) current 
sheets. In fact, such a state seems inevitable if the photospheric motion is large, so that 
the dislocation of the spines is large. It is evident from Figure 2 that reconnection shifts 
the inner spine to the left by transferring closed flux from overlying the left side of the 
polarity inversion line to the right. The amount of flux available for such transfer, 
however, is limited; consequently, if the ideal motions produce too large a dislocation of 
the spines, reconnection will not be able to realign them. 

Furthermore, there is no guarantee that reconnection will even preserve the basic 
spine-fan topology. 3D reconnection is likely to produce topologically complex structures 
so that the boundary between open and closed field becomes chaotic and the 
identification of a ID spine line or a 2D fan surface is no longer possible. This hypothesis 
seems even more likely if the closed bipole moves so that it encounters a large-scale 
coronal hole boundary. In that case the outer spine line would have to change from open 
to closed (or vice versa) and the fan would interact with the hole boundary. In order to 
determine the evolutionary topology and dynamics of 3D interchange reconnection, we 
calculate numerically two simple, but highly illustrative cases. In the first case (“open-to- 
closed”) an embedded bipole moves through an open field region and across a coronal 
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hole boundary, into a closed field region. In the second case (“closed-to-open”), we 
consider the reverse situation where a bipole moves from the closed field into the open. 
In Section 3 below, we describe the numerical model - initial conditions, driving flow 
field, and numerical grid. Section 4 discusses in detail the simulation results for the two 
cases. 

3. The Numerical Model 

We solve the standard set of 3D compressible, ideal MHD equations in spherical 
coordinates listed below, using the Adaptively Refined MHD Solver (ARMS) code 
(Welsch et al. 2005; DeVore, & Antiochos 2008; Lynch et al. 2008; Lynch et al. 2009; 
Pariat et al 2009): 

^ + V«pv = 0 
dt 

+ V • pvv + VP = — — (V xB \<B + pg 
dt 4ttV / 

on ^ 

+ W •Uv + P V »v = 0 

dt 

— = Vx(vxB\ 

dt V ' 


where all variables have their usual meanings. The internal energy density U is given by 
U = P/(y-l). The ratio of specific heats y is taken to be 5/3. The ideal gas law P = 2 ( p / 
m) k T is used as the plasma equation of state, where k is the Boltzmann constant and m 
is the proton mass. Gravity, given by g = - G Ms r /r , is included in the calculations, but 
its effects are small with regard to the interchange reconnection dynamics. The primary 
reason for adding gravity is to keep the plasma beta from becoming too large at large 
heights. 

The simulation domain consists of the spherical volume bounded below by the 
“photosphere” at r = 1 R s and bounded above by the source surface, which is taken to be 
at r = 3 R s . Within this domain, the initial magnetic field configuration is given by the 
analytic expression defined in section 2.1. The origin-dipole strength IM 0 I = 10 G, which 
yields a field strength of approximately 5 G at the photosphere far from the embedded 
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bipole. A single dipole with magnitude IM;I = 50 G, is placed below the surface at Ir) = 
0.9 Rs, the angular position of which varies between the two cases, although near the 
global coronal hole in both cases. Figure 3 shows the field for the case of the bipole 
initially in the coronal hole. 

A minor point to note is that we set the global dipole to be aligned with the y axis 
(0 = n/2, (|) = n/2) of the coordinate system rather than the vertical, as is the usual case. 
This implies that the coronal holes now occur centered about two points on the equator of 
our spherical coordinate system (at (j> = +/- n/2) rather than the coordinate poles (Figure 
3). Furthermore, we select the parameters of the embedded dipole so that it is located at 
and oriented along the coordinate equator, and impose photospheric flows that move the 
resulting embedded bipolar region along this equator toward a coronal hole boundary. 

The reason for this choice of geometry is that the poles have metric singularities in 
spherical coordinates, making them difficult to treat numerically, especially in 3D. The 
simplest and most effective procedure for dealing with these singularities is to remove 
from the computation domain a small conical region centered about each pole, (0 < 

1 1.25° in the north and 0 > 168.75° in the south), visible in Figure 3 as the “holes” in the 
Sun. We chose the magnetic and velocity fields so that all the structure and dynamics 
occurs at the equator, as far from these conical regions as possible. Note that there is no 
solar rotation in our simulation; hence, our choice of parameters for the magnetic field 
and flow fields corresponds only to a trivial rotation of coordinate axis and has no 
physical consequences. 

Since the initial magnetic field is potential, we set the initial plasma distribution to 
be spherically symmetric and in hydrostatic equilibrium: 



where the exponent u = Ro/H 0 = 1 1.66. The pressure scale height H 0 = 2 k T 0 / ( m go ). 

2 

The surface parameters are initialized to: To = 1 MK, P 0 = 1 dyne cm'". These plasma 
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profiles and parameters were selected so that the plasma (3 would be small throughout the 
domain. We find that [3 reaches a minimum of ~ 0.0325 inside the strong field of the 
bipole and an average value of less than ~ 0. 1 near the source surface; consequently the 
system is low-beta, as in the true corona. Furthermore, the gravitational energy of the 
plasma is small compared to the magnetic field energy. We emphasize, however, that 
although the system as a whole is low (3, the plasma pressure does play an important role 
in the evolution. Near the coronal null the plasma pressure dominates; therefore, the 
formation of the current sheets and the subsequent reconnection dynamics are critically 
dependent on the plasma evolution. 

Similar to the plasma beta, the Alfven speed varies considerably over the domain, 
but an average global Alfven speed can be defined as, 

y 1 2 Total Magnetic Energy 
1 y Total Coronal Mass 

With this definition, the Alfven speed in both simulations is approximately 400 x 10 5 
cm/s. An Alfven time of a little less than 1 15 minutes (ta ~ 6900 s) is similarly defined 
using a global length scale of 4 R s (about the length of the largest loops). 

At the lower boundary, the photosphere, we impose line-tied, no-flow-through (V r 
= 0) conditions. In both simulations, the embedded bipole is driven toward the coronal 
hole boundary by an incompressible surface flow applied at the photosphere, Figure 4. 
The flow field is constructed as a lst-order Fourier trigonometric series in the spherical 
angular coordinates. The azimuthal component, V^, is assumed to have cosine profiles in 
both colatitude (9) and longitude (<])) angular coordinates, and corresponding wave 
numbers that yield adjoining vortices (ke =1.0 and k,|, = 0.5). The polar flow component, 
Ve, is then calculated by applying the vanishing divergence condition for this 2D flow 
field. 
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where, k t = 0.5, and t H = 1.5 x 10 4 s = 2.17 ta. The magnitudes of these angular velocity 
components are set to be approximately an order of magnitude smaller than the average 
Alfven speed defined above; IV 0 I = 1.875 x 10 6 cm/s = 0.047 V A , and IV ( |,I = 5 x 10 6 cm/s 
= 0.125 V A . Note, however, that the driving speeds above are much smaller than the 
Alfven speed in the embedded bipole region, which is at least an order of magnitude 
larger than V A . In order to minimize transient wave effects as the motions start, the 
velocity magnitude has a shifted cosine profile in time. The flow is chosen to have a 
broad latitudinal range (0 H = 0.9 ti, 0 c = 0.5tc, Ql = O.I71) in order to minimize the 
distortion of the flux distribution within the embedded bipole as it moves across the 
photosphere (Figure 4). 

We use the velocity expressions above to describe the flows for both the case with 
the bipole initially in the coronal hole and the case with it initially in the closed field, 
except for a change in the longitudinal extent of the motions (and the obvious change in 
sign). In the first case, initially in the coronal hole, we set (j) H = 0.471, §c = 0.2tt, (])l = 0.0; 
whereas for the second case, we set (j) H = 0.7571, 4>c = 0.37571, 4>l = 0.0. These values for 
the flow parameters were selected so that the bipole would definitely cross the coronal 
hole boundary in both cases. 

At the top boundary, the source surface, we impose no-flow-through, free-slip 
conditions. The free slip conditions allow us to capture the physical distinction between 
open and closed field without having to incorporate in the simulations the added 
complexity of a solar wind. The key physical difference between closed and open field 
lines on the Sun is that the closed lines can contain long-lived stress (electric currents) 
whereas the open lines must be stress-free on non-transient time scales. Since closed lines 
have both ends line-tied at the photosphere, they can exert a finite stress at both ends. 
Open field lines, on the other hand, have only one end line-tied while the other is at 
infinity, so that any stress injected by photospheric motions will propagate out to infinity. 
(Note that the other end does not have to be actually at infinity, but only out beyond the 
Alfven point, which is typically at 20 solar radii, or so). We can capture this physical 
difference between open and closed by simply imposing free-slip conditions at the outer 
boundary. Field lines with both ends at the inner boundary are effectively closed, because 
they can hold stress indefinitely (time scales up to the dissipation time). Field lines with 
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one end at the inner boundary and one end at the outer are effectively open, because the 
free-slip condition implies that any stress injected on these lines will propagate out of the 
system on Alfvenic time scales. In addition to the free-slip, we impose a no-flow-through 
condition at the outer boundary so that we can preserve the open or closed property of a 
field- line under an ideal evolution. In our simulations, a field line can change from being 
open to closed or vice versa only as a result of reconnection and not by merely 
rising/sinking through the outer boundary. 

Finally, Figure 5 shows the numerical grid that is used for the simulations. We 
start with a base level consisting of 2 x 3 x 6 blocks distributed uniformly in r, 0, (j), with 
8 grid points per block. The initial minimum resolution is 3 levels refined above the 2 x 
3x6 base, and the grid is then refined a maximum of 6 levels over a volume 
encompassing the entire photospheric flow field and to a radius above the magnetic null 
point. The resolution at this highest refinement level corresponds to approx. 2.2 x 10 9 cm 
by 1.875° by 1.875°, which is much smaller than the scale of the embedded bipole or the 
flow field. Note that the grids are nearly identical for the two cases, except for minor 
adjustments due to the different initial position of the bipole and the latitudinal extent of 
the flow fields. In order to quantitatively compare the results of the two cases we have 
kept the grid fixed throughout the two simulations. 

4. Results 

4.1 Case 1: Bipole Convection from Open to Closed Field Regions 

Figure 3 shows the initial configuration for this simulation. The near-surface dipole is 
located at a latitude of 36.4°, which places the outer spine inside the coronal hole, but 
very near the coronal hole boundary (Figure 3). We chose this initial location so that the 
interaction between the embedded bipole field and the coronal hole boundary would 
occur before extreme distortion of the closed bipole field. The evolution for the 
convection of the bipole from the open to closed field regions can be considered to 
consist of 4 phases: 
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Phase 1. From t = 0 to t ~ 5900 s, the bipole moves toward the coronal hole boundary 
with evidence for only minor reconnection. Due to the finite grid of the simulation, some 
numerical resistivity is always present; therefore, if one examines field lines on a fine 
enough scale (less than the grid size), it is always possible to find some systematic flux 
transfer indicative of reconnection. The null point, however, remains almost undistorted 
during phase 1, and only weak currents (scale size substantially larger than the grid size) 
form there, so any reconnection is slow. The distance traveled by the inner spine during 

o 

this phase is ~ 34 x 10 cm, which is a small fraction of the scale of the bipole (the 

o 

diameter of the polarity inversion line in the direction of the motion is ~ 170x 10 cm, 
see Figure 3). As a result of the photospheric motions, the closed field region in front of 
the bipole is compressed, generating stresses on the open field. These stresses displace 
the inner and outer spines, exactly as in Figure 2, resulting in the eventual formation of a 
current sheet at the deformed null. 

Figure 6 presents a close-up of the null region, at t = 0, 5880, 7480, and 10000 
seconds. These are individual frames from a movie that is available in the electronic 
version of the paper. The white lines indicate initially open field lines and the yellow 
closed. Plotted on the vertical symmetry plane that bisects the bipole are filled contours 
of current density and 10 black contours of (3 with magnitude ranging from 1 to 100. This 
high (1 region corresponds physically to the null volume where the field is susceptible to 
strong distortion. It is evident from Figure 6 that the deformation of the null region stays 
small up through t = 5800, because the (1 contours remain approximately circular. The 
currents clearly build up as the bipoles motion progresses, but at this time they are still 
small compared to the currents produced by the driving motions. 

Phase 2. From t ~ 5900 s to t ~ 7500 s the continued motion of the bipole results in 
sufficient deformation of the null region that the structure of the currents there decreases 
down to the grid scale, and rapid reconnection occurs. This “interchange” reconnection 
exchanges the closed field of the bipole with the open field between it and the coronal 
hole boundary. It is most obvious in the movie, but it can also be seen in Figure 6. Note 
that the panel corresponding t = 7480 has substantially fewer white field lines to the right 
of the closed fan surface. 
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We find that once interchange reconnection turns on, it stays on and smoothly moves 
the outer spine through the open field and closer to the coronal hole boundary. There is 
little evidence for explosive dynamics such as bursty reconnection or large mass 
outflows. The dynamics produced by the interchange reconnection in this evolution are 
dramatically different than those in our simulations of breakout CMEs (e.g., Lynch et al. 
2008) or of coronal jets driven by magnetic twist (Pariat et al. 2009). The reason for this 
difference is that in the case of the CME and jet calculations, the photospheric motions 
are chosen so that that the magnetic stress is kept away from any separatrix surface. As a 
result, substantial free magnetic energy builds up inside the closed field volume until it is 
released by an explosive burst of reconnection, usually accompanied by some ideal 
instability or loss-of-equilibrium. In contrast, the large-scale translational motions of the 
simulation in this paper tend to move the bipole bodily, producing little magnetic stress 
inside its closed field. We find that only weak volumetric currents appear inside the fan 
and very little free energy is stored there. 

The motions do produce significant stress, however, on the large-scale field where the 
connectivity is discontinuous, the outer fan separatrix and outer spine. This stress leads to 
the formation of current sheets at the fan and null region, which are quickly dissipated by 
reconnection without large energy release or strong impulsive behavior, at least, for the 
magnetic Reynolds number of this simulation ~ 1000. Our result indicates that in order to 
obtain the large energy release to explain jets or plumes, for example, the closed field 
inside the fan would have to be stressed by small-scale photospheric motions as in Pariat 
et al. (2009) or emerge through the photosphere containing large stress. Both effects are 
almost certain to be true in the Sun due to the presence of subsurface convective flows 
and the photospheric granule and supergranule motions. 

Phase 3. Interchange reconnection continues until eventually the outer spine reaches the 
coronal hole boundary. At some instant around t ~ 7480 s the null of the closed field 
bipole lies exactly on the separatrix surface between open and closed field and, hence, the 
outer spine becomes a separator line that connects the bipole null and the null at the 
source surface. Of course, this is a singular event. At this time the coronal hole boundary 
can be taken to jump discontinuously from lying in front of the bipole to behind, so that 
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the fan bipole enters the main closed field region (Figure 6). Note that we see no 
evidence for any special dynamics during this period. The transition from the bipole 
being surrounded by open field to closed appears smooth. This result is to be expected, 
because the bipole field has such small scale that its interaction with the large-scale 
closed field just outside the coronal hole boundary is essentially identical to that of the 
open field inside that boundary. As far as the magnetic field of the bipole is concerned, 
there is negligible difference between the open and closed field regions. Furthermore, this 
result agrees with observations, which indicate that, in general, no special dynamics are 
seen at coronal hole boundaries (Kahler & Hudson 2002). 

Phase 4. During the final phase of the evolution, from t ~ 7500 s to t = 10000 s the bipole 
field moves steadily through the closed field by reconnecting with this flux. Note that 
although the total duration of the imposed flows is 15000 s, we end the simulation at t = 
10000 s; consequently the bipole is still being driven even at the end of the final phase. 
The reconnection during this phase is no longer of the interchange type, because it 
involves two closed field systems, but there appears to be little change in the dynamics. 
The current sheet at the deformed null region keeps increasing in length while decreasing 
in width (Figure 6), and the reconnection remains smooth with no apparent burstiness. 

We expect that if the bipole driving were to stop, the current sheet would decrease in 
length and the reconnection would eventually end, albeit with some residual currents left 
in the system. 

A critical issue is the topology of the open-closed boundary throughout this four-phase 
evolution. The quasi-steady models require that the reconnection maintains a smooth 
topology with well-separated open and closed field regions (Antiochos et al 2007). In 
order to determine the topology we have traced a dense sample of field lines from the 
source surface down to the photosphere and plotted their location there. Figure 7 shows 
the results for the open to closed simulation at three times during the simulation, and a 
movie showing the evolution at intermediate times is available in the electronic version 
of this paper. 
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The black region in each panel is the area on the photosphere that is magnetically 
connected to the source surface, in other words, the open field region. Also shown are 
the polarity inversion lines on the photosphere (thin black lines) and filled contours of B r 
at the photosphere, with red indicating strong negative and blue strong positive field. We 
note that in the first panel, at t = 5069 s, the bipole is completely surrounded by open 
field, so it is still in the coronal hole. The coronal hole forms an open corridor that 
extends around the negative polarity spot, but this corridor is fully connected at both ends 
to the main coronal hole open field region. This result shows that the mere observation 
of open field in strong active region magnetic fields does not constitute evidence for the 
validity of the interchange model. The quasi-steady models can easily account for such 
observations. 

As the bipole moves toward the closed field region, this open-field arch decreases 
in width due to interchange reconnection until by t = 7083, only a very thin corridor of 
open field remains (Figure 7). (The evolution is most clearly seen in the accompanying 
movie.) The key question is whether this corridor continues to be well-connected to the 
main open field region or whether it breaks up into disconnected segments. It does appear 
from Figure 7 that the corridor has breaks, but this is an artifact produced by the finite 
resolution of the numerical grid and the geometry of the photospheric flux distribution. 
Since the negative flux is concentrated into a spot just above the strong positive flux, it is 
relatively easy to find field lines that connect to this negative spot. This result is also 
evident in the initial potential field (Figure 1). If one draws a line connecting the centers 
of the positive and negative spots, the fan surface has a high density of field lines in that 
direction but low density in the perpendicular direction, so that the fan surface appears to 
have gaps in this perpendicular direction. We know from the analytic expressions, 
however, that the fan forms a smooth continuous surface. In topological terms, the 
reason for these apparent gaps is that the eigenvalues of the field Jacobian evaluated at 
the initial null point are highly asymmetric, so that the one corresponding to the 
eigenvector parallel to the center-to-center line is substantially larger than the eigenvalue 
for the perpendicular direction (e.g., Lau and Finn 1990). This asymmetry is maintained 
as the spots move and, hence, the open field corridor that eventually develops also 
appears to have gaps. However, when we plot field lines from the photosphere upwards 
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with very high resolution, we find that at t= 7083 there are always open field lines 
separating the closed flux that connects to the negative spot from the closed flux that 
connects across the equatorial inversion line. 

As the bipole moves toward the closed field region, the open field corridor 
continues to thin until eventually the outer spine coincides with the open-closed field 
boundary boundary, so that the corridor achieves singular width. Since our simulation has 
finite temporal and spatial resolution, we cannot capture this critical event when the 
corridor is singular. It is possible that near this time the open field corridor breaks up into 
discontinuous pieces, because the deformation of the null and the presence of current 
sheets there cause the outer spine to deform to a sheet-like structure and the fan to some 
fractal volume. We do not see such topologies in the simulation, the outer spine remains 
ray-like, but this may be due only to the finite resistivity inherent to our numerical code. 
Even if such singular topologies do occur, we expect that their structure would be only of 
order the dissipation scale and, consequently, disappear quickly. Our simulation shows 
only a smooth topological transition for the bipole as it moves from the open to closed 
regions, in good agreement with the results of the quasi-steady models (Figure 7 and 
accompanying movie). 

4.2 Case 2: Bipole Convection from Closed to Open Field Regions 

The closed-to-open case is, for the most part, closely analogous to the open-to-closed 
evolution. The near-surface dipole is initially located at the latitude of 36.0° placing the 
outer spine inside the closed field, very near the coronal hole boundary in order to 
minimize distortion during the interaction (Figure 8a). Again, we organize the evolution 
of the bipole from the closed to open field regions into 4 phases: 

Phase 1. From t = 0 to t ~ 3758 s, the bipole moves toward the coronal hole boundary 
with little reconnection or current sheet formation. The distance traveled by the inner 

o 

spine during this phase is ~ 80 x 10 cm, about half of the dipole polarity inversion line 
diameter. The photospheric motions expand the entire global closed field region, 
generating magnetic field stresses behind the dipole. The inner and outer spines separate 
as a result of these stresses, eventually deforming the null and generating a current sheet. 
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Figure 8 shows the evolution similar to the open-to-closed case, at t = 0, 3758, 9273, and 
10000 seconds - movie available in the electronic version of the paper. Clearly, from 
Figure 8, the deformation of the null region stays small up through t ~ 3758 s as the (3 
contours are still approximately circular. The currents within the null region build up as 
the bipole motion progresses, but they are still small compared to the driving motion 
currents. 

Phase 2. From t~3758 stot~ 9273 s, the continued motion deforms the null region, 
decreases the current structure to the grid scale, and initiates rapid reconnection. Though 
not strictly interchange reconnection because the bipole is embedded in a globally closed 
field region, reconnection interchanges the closed flux inside the bipole fan separatrix 
with the large-scale closed field. Once again we find that the system evolves by 
continuous reconnection, smoothly shifting the outer spine through the embedding field, 
with little evidence of bursty reconnection or large material outflows. As above, only 
weak volumetric currents appear inside the fan and very little free energy is stored there. 

An important difference between this case and the open-to-closed case is that the 
displacement required for the bipole to cross the coronal hole boundary is much larger 

o 

than before. The required displacement is ~ 680 x 10 cm, nearly 4 times the bipole 
polarity inversion line diameter. This result is due to the difference between the response 
of open field and closed field to photospheric stressing. Since the open field is free to slip 
at the source surface, significant compression stresses do not build up between the front 
of the bipole and the coronal hole boundary. For this case, the photospheric motions 
stress primarily the closed field, which is line-tied at both footpoints. Consequently, the 
stress at the null and fan surface originates from behind the bipole as a result of the 
stretching of the closed field there. However, the eventual results of this stress are the 
same: current sheets form along the fan and deformed null region and dissipate quickly 
by reconnection without large energy release or strong impulsive behavior. 

Phase 3. At some time around t ~ 9273 s, reconnection between the bipole flux inside 
the fan surface and the external field shifts the outer spine line to the coronal hole 
boundary, so that the boundary jumps discontinuously across the bipole fan surface 
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(Figure 8). Again, this singular topological transition appears smooth, showing no 
evidence of any special dynamics. 

Phase 4. During the final phase of the evolution, from t ~ 9273 s to t = 10000 s the bipole 
field moves steadily through the coronal hole by reconnecting with the open field. The 
reconnection during this phase is true interchange, because the bipole is now embedded 
in the open field region. The current sheet aspect ratio continues to increase at the 
deformed null region (Figure 8), and the reconnection remains smooth. We expect that if 
the bipole driving were to stop, reconnection would eventually dissipate the current sheet. 
Since the motion is now within the open field, any helicity injected by the photospheric 
motions may escape the system allowing a realignment of the inner and outer spines. For 
this case it is possible that the system can achieve a true minimum-energy, potential state, 
except perhaps for any volumetric currents deep inside the closed bipole field. 

Finally, we find that the evolution of the magnetic topology (Figure 9), is essentially 
identical to that above. Initially, the parasitic spot is completely surrounded by closed 
flux. At some point near t = 8902 s, the bipole is so close to the coronal hole boundary 
that the outer spine shifts its global topology, and a very thin open field corridor forms. 
Using the same arguments as in the open-to-closed case, the open field corridor is 
expected to be well connected to the main coronal hole even though in the figure it 
appears to have breaks. Once the motion is completely inside the open field region, the 
corridor continues to thicken as flux is transferred across the bipole fan surface (Figure 
9). The sequence shown in Figure 9, therefore, is simply the reverse of that in Figure 7. 

5. Conclusions 

The results described above have a number of implications for understanding the 
corona/heliosphere and interpreting observations. The first and, perhaps, most important 
conclusion is that the basic topology of the interchange process is that of the closed field 
of a bipole interacting with surrounding open field, as in Figure 2. The reconnection 
occurs at the fan surface, primarily at the null. Note that the topology is continuous and, 
hence, it is not valid physically to assume a picture in which reconnection takes place 
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between an isolated open and closed field line. The difference between the continuous 
topology of Figure 2 and the often-used discontinuous picture may seem minor, because 
in both models the open field undergoes a jump in footpoint position as a result of 
reconnection. The key point, however, is that in the continuous model the reconnection 
releases energy only after a large current sheet forms. If the reconnection at the null is 
highly efficient, the open field will smoothly transfer from one side of the bipole to the 
other with no heating or mass acceleration. We conclude that even though copious 
interchange reconnection may be present in coronal holes, this reconnection may not play 
a significant role in the heating and acceleration of the wind. A robust treatment of the 
reconnection, especially of the effective resistivity, is required in order to evaluate the 
importance of the interchange process to the formation of the solar wind. 

From our simulations, the magnetic topology appears to remain fairly smooth 
throughout the interchange reconnection process, even when a bipole crosses a helmet 
streamer boundary. Our results, therefore, constitute strong support for the quasi-steady 
models in which the large-scale field evolution can be approximated as a sequence of 
topologically smooth quasi-steady states. One aspect of this general result is that the 
uniqueness conjecture (Antiochos et al. 2007) appears to hold even during interchange 
reconnection. We see no evidence for disconnected coronal holes as the bipole evolves, 
which argues against the basic assumptions of the interchange model. It should be 
pointed out, however, that we have calculated the evolution of only a single bipole 
moving in a simple trajectory. The interchange model of Fisk et al (1999) inherently 
assumes that the evolution of the open field is dominated by its reconnections with a 
dynamic complex of random bipoles, as in the magnetic carpet. It may well be that if a 
sufficiently complex distribution of bipoles is present, then key features of the 
interchange model, such as open-flux diffusion, start to become valid. 

Another important result from the simulations above is that unlike the source- 
surface model, the field does not remain current-free during the evolution. Large currents 
do form in the corona in response to the photospheric motions, and these currents are 
long-lived. Consequently, the position and geometry of the open-closed boundary will be 
different than that calculated from the source surface model, which is important for 
comparison with observations, but the model should still be qualitatively correct. Given 
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the proper boundary conditions, the MHD models can, in principle, calculate the field and 
currents precisely, but determining such boundary conditions from available observations 
may not be possible. 

From the viewpoint of comparison with observations, a key result from our work is 
that the reconnection does not produce bursty dynamics. Energy is released during the 
reconnection, primarily as mass flows, but the release does not show the type of 
impulsive behavior that we have seen in previous simulations of coronal reconnection 
(e.g., Karpen et al. 1998; Pariat et al 2009). There are several possible reasons for this 
difference. First, the photospheric driving in the present simulations does not impart 
significant shear or twist to the field, unlike the case of our coronal jet model (Pariat et al 
2009). Consequently, the energy released by the interchange reconnection is only the 
energy in the local current sheet at the null region, not the free energy of currents in a 
large volume. Second, the reconnection is fully 3D, which makes it more difficult for 
magnetic islands to form. In our 2.5D simulations (Karpen et al. 1998) we found that 
much of the burstiness was due to the random formation and expulsion of magnetic 
islands in the reconnecting current sheet. Another problem is that, even with adaptive 
mesh refinement, the 3D calculations simply have less numerical resolution that the 2.5D 
runs, so that the effective diffusion is larger. It may well be that with sufficient resolution, 
the 3D runs would produce current sheets with very large aspect ratios and the 
reconnection dynamics would more resemble those of the 2.5D results. 

On the other hand, observations of bipoles in coronal holes do not appear to show 
strong continuous dynamics; consequently, our results above may hold, in general. 
Coronal jets are relatively rare events. In most cases embedded bipoles are associated 
with long-lived plumes, which require a quasi-steady heating (e.g., Wang and Muglach 
2008). The type of reconnection that we obtain in the simulations above would be 
compatible with this type of heating. It is difficult, however, to compare our results 
quantitatively with plume observations, because we do not include a proper treatment of 
the plasma energetics. Our initial conditions assume a spherically symmetric density 
distribution, but on the Sun, the densities inside the closed field region of a coronal-hole 
bipole are clearly much larger than the density of the surrounding open field. Interchange 
reconnection between the closed and open field will release this high-density (and 
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pressure) material onto open field lines, which would then expand outward and, perhaps, 
form a plume. Further calculations with more realistic plasma energetics, including 
radiation and thermal conduction are needed in order to test such a model. 

It is tempting to conjecture that this process of releasing the closed-field plasma of 
embedded bipoles onto open field is the origin of the slow wind. The problem, however, 
is that closed bipoles are observed to occur throughout coronal holes, so that if this 
process were important, slow wind would be observed far from the heliospheric current 
sheet. The fact that slow wind is seen only near the current sheet, within 15° or so, 
implies that the origin of the slow wind must be associated with coronal hole boundaries. 
The results shown in Figures 7 and 9 indicate that if there are many bipoles moving 
randomly in response to photospheric motions, the coronal hole boundary is likely to 
consist of a complex dynamic web of open-field corridors. Such a dynamic topology 
would blur the distinction between the interchange and quasi-steady models, at least in 
the vicinity of coronal hole boundaries. We conjecture that a mixing of the two types of 
models is, in fact, the key to understanding the origins of the slow wind. 

Acknowledgements 

This work was supported, in part, by the NASA HTP, TR&T, and SR&T Programs. All 
the high performance computing capabilities were provided by the DoD HPCMP. J. K. 
Edmondson gratefully acknowledges support of a NASA GSRP grant for his PhD 
research. 


- 26 - 




Figure 1. Magnetic topology of a small near-surface dipole and global dipole, a) Colored 
contours show magnitude of radial field at photosphere, the two white curves indicate 
polarity inversion lines (radial field vanishes). The yellow field lines above the surface 
correspond to streamer belt closed flux and the white field lines to the open, coronal hole 
flux that maps to the source surface, b) Close-up of the field near the embedded bipole 
showing the outer fan field lines and spine. 
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Figure 2. Interchange reconnection schematic: A) Initial field configuration. B) Stressed 
field configuration. C) Current sheet & reconnection jets D) Interchange reconnection 
flux exchange 
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Figure 3. Global magnetic field configuration. Open field coronal hole regions are shown 
in white. The closed field, streamer belt region is shown in yellow. The spine fan 
topological characteristics of the embedded bipole are shown in red. 



Figure 4. Surface velocity field vectors. Color scale: red indicates zero velocity 
magnitude, purple indicates spatial extent of flow field. 
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Figure 5. Numerical grid. Top panels: Grid refinement in the radial direction. Bottom 
panels: Grid refinement across the surface. Note, the initial minimum refinement is 
refined to 3 levels above the base 2x3x6 blocks. 
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Figure 6. Open-to-closed evolution, a) t = 0 s; initial configuration, b) t = 5,880 s; current 
sheet formation, c) t = 7,480 s; global topology change of external spine, d) t = 10,000 s; 
final configuration. 
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Figure 7. Open flux mapping evolution on photosphere. Open-to-Closed case at t = 
{5069, 7083, 8960} seconds. 
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Figure 8. Closed-to-open evolution, a) t = 0 s; initial configuration, b) t = 3,758 s; current 
sheet formation, c) t = 9,273 s; global topology change of external spine, d) t = 10,000 s; 
final configuration. 
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Figure 9. Open flux mapping evolution on photosphere. Closed-to-Open case at t = 
{7752, 8902, 9915} seconds. 
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